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Abstract 

We use gauge theory/string theory correspondence to study the bulk viscosity of 
strongly coupled, mass deformed SU{N^ A/" = 4 supersymmetric Yang-Mills plasma, 
also known as A/" = 2* gauge theory. For a wide range of masses we confirm the bulk 
viscosity bound proposed in [1]. For a certain choice of masses, the theory undergoes 
a phase transition with divergent specific heat cy ~ |1 — Tc/T|~^/^. We show that, 
although bulk viscosity rapidly grows as T — T^, it remains finite in the vicinity of the 
critical point. 
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1 Introduction 



In [2] Maldacena proposed that SU{Nc) Af = A superconformal Yang-Mills (SYM) 
theory on R^'^ is dual to type IIB string theory on (a Poincare patch of) AdS^ x 
with Nc units of the Ramond-Ramond five-form flux through the S^. Assuming that 
such a duality holds exactly at a superconformal flxed point, it must also hold for 
any relevant (in the infrared) deformation of a fixed point. Specifically, we might 
generate non-conformal examples of gauge theory/string theory correspondence by 
simply mapping the mass deformation on the SYM side to the string theory side. 
For infinitesimal supersymmetric mass deformations this was done in [3] using the 
operator /state correspondence in AdS/CFT [4, 5]. Extending an infinitesimal mass 
deformation to a finite one proved to be extremely difficult. In fact, only one particular 
deformation was shown to be integrable in the 't Hooft limit Nc — > oo, Qym (with 
A = NcQym kept constant), and for large 't Hooft coupling A ^ 1 [6] (PW). In [7] 
it was shown^ that PW massive deformation is dual to giving the same mass to two 
chiral multiplcts of the parent J\f = A SYM at a specific point on its Coulomb branch^. 
Such a mass deformation breaks supcrsymmetry down to A/" = 2. The resulting gauge 
theory, commonly referred to as A/" = 2* gauge theory, can be solved nonpcrturbatively 
[9]. The agreement [7] between the dual gravitational description of the massive M = 2 
gauge theory [6] and its exact field-theoretic solution [9] provides a highly nontrivial 
check of gauge theory/string theory correspondence in a non-conformal setting. 

Once the duality is established for a supersymmetric ground state, following [10] it 
can be extended to correspondence involving a thermal equilibrium state of the gauge 
theory: a thermally equilibrium state of a gauge theory is dual to a Schwarzschild black 
brane solution in the supergravity description, with the temperature T given by the 
Hawking temperature of the black brane. In the context of A/" = 2* theory there is an 
interesting subtlety: the supersymmetric mass deformation involves deformation of the 
parent SYM by operators of different canonical dimensions, a dimension-two operator 
for the bosonic components of massive chiral multiplcts and a dimension-three operator 
for the fermionic components of massive chiral multiplcts. Such operators are mapped 
to different scalar gravitational modes of the effective five- dimensional gravitational 
description [6]. The coefficients of the non-normalizable modes of these scalars encode 

-'^ Related discussion appeared in [8]. 

^Extending the mass deformation duality to all Coulomb branch of the M = A SYM is a difficult 
open problem. 
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the bosonic nib and fermionic ruf masses respectively [3]. Although the vacuum state 
supersymmetry requires nib = irif, a, thermal state breaks the supersymmetry anyway; 
thus, we can study a phase diagram of A/" = 2* gauge theory with rub 7^ [H]- A 
full ten-dimensional type IIB supergravity solution dual to A/" = 2* gauge theory for 
generic {T.nib.nif) was constructed in [11]. It was shown there (see also [12]) that 
any such state (including small fluctuations about it) can be described within effective 
flve-dimensional gauged supergravity presented in [6]. 

Supergravity equations of motion derived in [11] describing a thermal state of 
J\f — 2* gauge theory are too difficult to solve analytically. If fact, in [11] these equa- 
tions were solved only in high-temperature limit m;,/T -C 1 and nif/T^ 1. Additional 
analysis of the M — 2* phase diagram required numerical work. For two special cases, 
i.e., (susy) = {nib = ''^/j and (bosonic) = {nif = 0,mft 7^ 0}, such numerical analy- 
sis were performed in [13]. Using the holographic renormalization of the theory [14], 
the energy density S, the free energy density and the entropy density s of A/" = 2* 
strongly coupled plasma was computed. It was shown in [13] that the basic thermody- 
namic relation = S — Ts is satisfied exactly, while the first law of thermodynamics 
d8 = Tds is exactly satisfied for the high temperature analytic solution of [11], and 
is satisfied numerically for [susy) and (bosonic) thermal states with an accuracy of 
~ 0.1% and ~ 0.01% correspondingly. The latter provides a highly nontrivial check 
on our identification of the bosonic and fermionic masses in dual supergravity (see [13] 
for details). 

An interesting critical point was found in the numerical analysis of the (bosonic) 
thermal state of the M = 2* plasma [13]: for T < ^ mb/2.29(9) M = 2* plasma 
becomes unstable with respect to energy density fiuctuations. Specifically, precisely at 
T — Tf. the speed of sound waves squared vanishes. A pcrturbative instability of this 
type is a defining feature of a second order phase transition. Since oc (T — Tc)^/^, 
the specific heat Cy diverges as [1 — Tc/T|^^/^, suggesting that such a critical point 
is in the universality class of the mean-field tricritical point. Physically, the existence 
of pcrturbative instability in = 2* plasma at low temperatures is not surprising. 
Indeed, once nib the supersymmetry is broken, and the theory is guaranteed to 

be stable only at high temperatures. It was conjectured in [13] that J\f — 2* plasma 
would have a critical point 




(1.1) 



as long as ly < 1. 
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Shear viscosity of A/" = 2* plasma was computed in [15, 14]; it was shown to satisfy 
the universal bound [15, 16, 17] 

s Air 

Transport properties ol M — 2* plasma were further studied in [18]: 

■ equations of motion describing sound quasinormal modes oi H — 2* black brane 
with arbitrary momentum q = |^/(27rT') and frequency tt) = ijj/{2t:T) were obtained; 

■ these equations were solved analj^ically in the hydrodynamic hmit q ^ 1 , tn <^ 1, 
and for high temperatures T ^ {m^, m/}; 

■ following [19], the dispersion relation for the lowest quasinormal mode of the J\f — 2* 
black brane was identified with the dispersion relation of the sound waves in strongly 
coupled J\f — 2* plasma: 



tt) 



= c.q-2.z^Q + |-) q^ + 0(q^), (1.3) 



where ( is the plasma bulk viscosity; 
■ to leading order in {ruf /T,mb/T} it was found that 

in precise agreement with the speed of sound obtained from the analytic equation of 
state oi J\f = 2* plasma at high temperature^ 

4--%, (1.6) 



using (1.2), (1.4), to leading order in (| — cfj it was found that 

'1 



where /3F ^ 8.001; 



mf=0 



16 V3 



37r/3y /I 



O 



(1.6) 



'/ mb=Q 

where ^ 0.66666 [20]. 



2 V3 



(1.7) 



^This provides a highly nontrivial consistency check on our analysis of the quasinormal modes. 
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In this paper we extend analysis of [18, 1] to general mass deformations, i.e., apart 
from cases of {hosonic) <^ i/ = and {susy) u = 1 thermal states oiN = 2* plasma, 
and for wide range of temperatures. Our goal is twofold: 

■ first, we would like to test the bulk viscosity bound conjecture of [1] in more general 
setting: 



for < 1/ < 1. 

The paper is organized as follows. In the next section we outline the equations of 
motion and the boundary conditions for the J\f — 2* black brane background and its 
hydrodynamic quasinormal mode. The results of our extensive numerical analysis are 
presented in section 3. We conclude in section 4. 

Most technical details are omitted due to their complexity. All equations of motion, 
their analytic asymptotic solutions, as well as raw numerical data is available from the 
authors upon request. We use numerical techniques developed in [21]. 

2 J\f = 2* black brane and its hydrodynamic quasinormal mode 

We closely follow [13] in discussion olM — 2* black brane background, and we (mostly) 
follow notations of [18] in discussion of its quasinormal modes. 

2.1 Effective action 

The effective action of the five-dimensional gauged supergravity describing J\f — 2* 
black brane thermodynamics/hydrodynamics is given by 




(1.8) 



■ second, we would like to compute 



(1.9) 



critical ^ T=T^{v) 




(2.1) 
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where the potential^ 



V 



1 

16 



1 fdW 
2)\ da 



(dW\ 



(2.2) 



(2.3) 



is a function of a and x, and is determined by the superpotential 

W = -e"^" - ^e^" cosh(2x) . 
In our conventions, the five-dimensional Newton's constant is 

2.2 M —2* black brane background 

We parameterize the background metric of the M — 2* black brane as 

dsl = C2{xf {-{I- xfdt^ + dxl + dxl + dxl) + Qxxix) dx^ , (2.5) 

where the radial coordinate x e [0, 1], with x — > 0+ corresponding to the AdS^, bound- 
ary and X — > 1_ corresponding to a regular Schwarzschild horizon of the black brane. 
The metric (2.5) is supported by nontrivial profiles of two scalar fields: 

Pq{x) = e^" , c{x) = cosh(2x) • (2.6) 

Notice that we redefined the gauged supergravity scalars as in (2.1) — this is done in 
order to speed up numerical integration. We further introduce 



C2 = e^, 



A = In ^3 - ^ ln(2a:; - x^) + a{x) , 



(2.7) 



where ^3 is related to the Hawking temperature of the black brane as follows [13] 

f 

T = lim e-^'*^^) . (2.8) 



The asymptotic solution of {pq, c, a} near the boundary, x — > 0+, takes the form 



[13] 



P6 



a/2 



(6pio + 6pii \nx) + O (a: In^ x) , 



c =1 + 12z/ x^/^ii + 24i/pii X (xio + i^Pii + 2upn lnx) + {x^/^ In^ x) , (2.9) 
2 

a = - -upux^/^ + 0{x In^ x) , 
o 



'*Wc set the five-dimensional gauged supergravity coupling to one. This corresponds to setting the 
radius L of the five-dimensional sphere in the undeformed metric to 2. 
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and near the horizon, y = 1 — x — > 0_|_ [13] 

p, = pl + 0{y^), c=4^ + (^(y'), a = ao + a,y^ + 0{y'). (2.10) 

In (2. 9), (2. 10) wc indicated expUcitly only terms necessary to unambiguously determine 
the asymptotic black brane geometry for a fixed set {n^v}. pu is related to p = ^ 
by [13] 

Notice that for a fixed set {//, u}, J\f — 2* black brane geometry is specified by six 
parameters: 

{^,1/} =^ |pio,Xio,Po,Xo,ao,ai| , (2.12) 

which is precisely the number of parameters needed to uniquely determine the solution 
of three second order equations of motion for {pg, c, a}. These parameters are functions 
of {p,iy}. 

We use numerical techniques developed in [21] to generate data sets 

^background — 

which can be used to study the thermodynamics oi Af — 2* black branes as detailed in 
[13]. The data sets Sbackground are available from the authors upon request. We obtain 
the following results. 

■ For u > 1, in agreement with the conjecture in [13], we did not find a critical 
point in A/" = 2* plasma down to ~ 10, which corresponds to temperatures of order 
T ~ mb/10. As ioT u — 1 in [13], the low temperature {p ~ 10) thermodynamics of 
J\f = 2* black brane can be well approximated by the following equation of state 

To^-NlT^e-^, (2.14) 

where rueff = meffiy) ~ m^/lO. 

■ For 0<z/<1,jV = 2* brane thermodynamics has a critical point — Tdv), such 
that 

- = = -i^ oc ±(r - Te) V2 , (T - Te) « % . (2.15) 
Cv oc 
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2.3 Quasinormal sound mode 

Sound quasinormal modes in the J\f — 2* black brane geometry involve coupled fluctu- 
ations of gauge-invariant metric fluctuations Zh and the gauge-invariant fluctuations 
of two scalar fields {Z^, Z^}, see [18] for details. The spectrum of quasinormal modes is 
determined [19] by imposing on {Zh, Z^, Z^} an incoming wave boundary condition at 
the horizon, and requiring vanishing of the non-normalizable modes for {Zh, Za, Z^} 
near the boundary. In the hydrodynamic limit tn — > 0, q — > with ^ kept fixed, this 
leads to the following perturbative expansions 

Zh =(1 - x)-^"' (^zh,o + ^q zh,i + 0{q' 

Z^ =(1 - (z^,o + ^(\ ^a,i + 0(q') ) , (2.16) 



(2.17) 



Z^ =(1 - x)-^^ \^z^,o + ^ z^,, + 0{q')j , 

with the following boundary conditions on {zH,i, Za,i, z^^i}: 

lim zho — i , lim zhi — , lim Za % = lim z^ i — finite , 

_ X— >1_ X— »1_ X— »1_ 

ZH,^ = 0{x) , Za,i = O {x^/^) , z^,, = O {x^/") , as x^0+. 
Additionally, we parameterize dispersion of the lowest quasinormal mode as 

to^^q(3,-'-q' (32 + 0(q'), (2.18) 

where (3i = f). In the conformal case, /i — * (with t/ = 0) we expect [22] 

lim /3i(//, 0) = 1 . (2.19) 

Identifying (1.3) with (2.18) and using the universal result for the shear viscosity 
(1.2) we find 



^-cij=-{l-/3t), ^=3(/52-l). (2.20) 
2.3.1 Leading order in the hydrodynamic approximation 

To leading order in the hydrodynamic approximation, wave functions of the gauge- 
invariant fluctuations {zh,o, Za,o, Zx,o} satisfy the following equations 

~^'h,0 + ^101 z'ha + C\02 Zh,0 + Cl03 Za,0 + Cl04 ^x,0 ) 

~Za,0 + ^201 z'^ Q + C202 z'h^o + ^203 ^Hfi + ^204 Zafl + ^205 %,0 j (2-21) 
=z'^^o + C301 z'^^Q + C302 ^jj^o + ^303 Zh,0 + C304 ^^a.O + ^305 ^x,0 ) 



where connection coefficients Cjoj ^-re nonlinear functionals of {pq, c, a} witli explicit 
dependence on x and /?i2 = Pf: 



{p6,c,a}; x; (3i2 



[2.22] 



Using (2.9), (2.10), for each set Sbackground we construct the asymptotic solution of 
(2.21). Near the boundary, x — > 0+ we find 

Zh,0 = X ^^ 2,0 + ^ (^^ 1^^) ' 

/ \ (2 23) 

Za,o = x^' 4°!,o + 0{x\nx), z^,o = x^' i^x'/' 4°i,o + 0{x\nx)j, 

and near the horizon, y = 1 — x ^ 0+ 

z„,o = 1 + 0{y') , Za,o = ?r + ^(y') , ^x,o = 4°^ + 0{y') . (2.24) 
Thus, altogether we have six new parameters: 

<j jo (0) (0) (0) (0) (0)1 /'O OCC^ 

'^background ^ 1 Pl2) -2^i/,2,0' '^a,l,0' -^XA-O' % '"^0 f' \^..iO) 

precisely what is necessary to construct a unique solution for {zh,q-i z^fi, z^^} for a 

given Siia^i^gyQund' 

We use numerical techniques developed in [21] to general data sets 

O =/c ■ R (0) (0) (0) (0) (0)1 ('9 Ofi^l 

'-'sound — '-'background) Hl2-i ^i7,2,0' ^a,l,0' ^x,l,0' % ' -^0 ( ' l^Z./OJ 

The data sets S sound a-re available from the authors upon request. 

2.3.2 The first subleading order in the hydrodynamic approximation 

To the first subleading order in the hydrodynamic approximation, wave functions of 
the gauge-invariant fiuctuations {zh,i, Za,i, satisfy the following equations 

=Zh i + Cm z'jj^i + Cii2 Zh,1 + C113 Za,l + Cii4 Z^^i + C115 2;^_o + ^116 Zh,0 

+ C211 ]^ + C212 z'jj j^ + C213 + C214 ^a,l + C215 ^x,l ^216 ^a,0 + ^217 ^iJ,o 

+ C2I8 ^if,0 + C219 Za,0 + C2110 ^x,0 ) 
=^^,1 + C311 z'^ i + C312 % 1 + C313 Zh,1 + Csu Za,l + C315 + C316 z'^^Q + C317 

+ C318 ^H,o + C319 2;q,o + C3110 ^x,o ) 

(2.27) 



where connection coefficients Cnj are nonlinear functionals of {pq, c, a} with exphcit 
dependence on x and {/9i2,/92}: 



{P6,c,a}; x; {f3i2, P2} 



(2.28) 



Using (2.9), (2.10), (2.23), (2.24), for each set Sgound we construct the asymptotic 
solution of (2.27). Near the boundary, a; — > 0+ we find 

zh,i = X z%]^ Q + O [x^ Inx) , 

/ \ (2 29) 

Za,i = x'^' 4^!,o + 0{x\nx), z^,, = x'/' (x'/' 4^1,0 + 0{x\nx)y 

and near the horizon, y = 1 — a; — > 0+ 

ZH,i = + 0{y') , = gj^) + 0{y') , z^,, = x« + O(y^) . (2.30) 

Thus, altogether we have six new parameters: 

fD M) Jl) Jl) Al) Jl)\ (n 0-1 N 

y'^T Zfj2fi, z^ i Q, z^ i Q,qQ ,Xq ?, l^'-J-LJ 



^sound 



precisely what is necessary to construct a unique solution for {zh,i, Za,i, z^^i} for a 
given S sound- 

We use numerical techniques developed in [21] to general data sets 

o — J Q . P W (1) (1) (1) (1)1 (0 "io^ 

<^ attenuation — S sounds P2, 2,0' ^a,l,0) ^x>l,0' ^0 ; -^0 f ■ [Z.OZj 



The data sets Sattenuaticm are available from the authors upon request. 
3 Bulk viscosity of jV = 2* plasma 

Given Sattenuation wc havc all the necessary data to study bulk viscosity of strongly 
coupled Af — 2* plasma. Of primary interest to us are the bulk viscosity bound 
conjecture [1], and the behaviour of bulk viscosity near the phase transition. 

3.1 Bulk viscosity bound 

Fig. 1 represents the ratio of bulk to shear viscosities (/r] as a function of (| — c^) for 
M = 2* plasma with mass deformation parameter v ranging from v = 0.2 to u = 0.9 
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0.02 0.04 0.06 0.08 0.10 0.12 0.14 ^3 " 

Figure 1: Ratio of viscosities - versus the speed of sound in A/" = 2* gauge theory plasma 

771 

with mass deformation parameter u = & [0.2,0.9], with intervals Au = 0.05. The 
solid line represents the bulk viscosity bound (1.8). 

with intervals of Au = 0.05. Different color sets of points represent different values 
of u. The solid blue line represents the bulk viscosity bound (1.8). We verified that 
in the vicinity of = | (which corresponds to a high-temperature regime of A/" = 2* 
plasma), the results are in excellent agreement with (1.6), (1.7). Although in this paper 
we truncated the plots to a near-conformal regime^, available data sets Sattenuation allow 
us to study the behaviour of bulk viscosity near the critical point Tc(z/). Much like for 
the 1/ = case discussed in [1], for each value of z/ < 1 we observe a rapid, power- law 
like, growth of C/r/ in the vicinity of the critical point. Nonetheless, this ratio is finite 
precisely at T = Tc (see Fig. 3 below). The results are qualitatively identical to the 
u = Q case discussed in [1]. 

Fig. 2 represents the ratio of bulk and shear viscosities C/r/ as a function of (| — c^) 

'^The rapid growth of bulk viscosity indicates that if any violation of the bound (1.8) would occur, 
it would occur in a near-conformal regime. 
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c 

V 




0.02 0.04 0.06 0.08^3" 

Figure 2: Ratio of viscosities - versus the speed of sound in A/" = 2* gauge theory plasma 
with mass deformation parameter i/ = G {1, 1.05, [1.1, 2.1]a;/=o.15 [3, 6]ai/=i}- The 

b 

sohd hne represents the bulk viscosity bound (1-8). 

for M = 2* plasma with mass deformation parameter v = {1, 1.05}, also in the range 
from z/ = 1.1 to z/ = 2.1 with intervals of Au = 0.1, and in the range from z/ = 3 
to z/ = 6 with intervals of Au = 1. Different color sets of points represent different 
values of z/. The solid blue line represents the bulk viscosity bound (1.8). Again, we 
verified that in the vicinity of = | (which corresponds to a high-temperature regime 
of A/" = 2* plasma), the results are in excellent agreement with (1.6), (1.7). Notice 
the accumulation of the points as data sets approach the bound. Much like for the 
z/ = 1 case discussed in [1], this regime corresponds to a low-temperature Af = 2* 
plasma regime. Specifically, there we have rrib/T ~ 5 • ■ ■ 10. Much like in [1] we can 
extrapolate the speed of sound and the viscosity ratio to T 0. We find that the 
endpoints of such extrapolations land on (or slightly above) the bulk viscosity bound 
line. 

It is very interesting to explore the thermodynamics of A^ = 2* plasma (for u > 1) at 
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6.8 



6.7 



6.6 



6.5 



0.2 



0.4 



or6 



0.8 



Figure 3: Ratio of viscosities ^ at T = Tc versus the mass deformation parameter u of 
M = 2* plasma. 



extremely low temperature^, ideally, as T ^ 0. So far we have no indication that such 
a zero temperature limit would be singular. On the other hand, we can not exclude 
the presence of some exotic instabilities/phase transitions at very low temperatures, 
see [11] for more details. 

3.2 Bulk viscosity at a critical point 

We already mentioned that M = 2* plasma with mass deformation parameter < 
u < 1 undergoes a phase transition, which appears to be in the universality class of 
the mean-field tri critical point. At such a phase transition, specific heat diverges as 
|1 — Tc/T|~^/^. Fig. 3 represents the ratio of bulk to shear viscosities C,/ri as a 



2 

m 



function of z/ = — j at the critical temperature T = Tc{h'). We took the z/ = result 



from [1]. Notice the rapid growth of - 



as V 1_. Actually, such a behavior is 

critical 



^Our current numerical algorithms become very slow as low temperatures. 
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expected if z/ = 1 is the end point at which transitions cease to occur, and for > 1 
the ratio oi (/rj is relatively close to the bound (1.8) down to rather low temperatures. 

4 Conclusion 

In this paper we study the bulk viscosity of A/" = 2* gauge theory plasma at strong 
coupling as a function of temperature and for various masses, u = mj/ml. In all 
cases, we find that the viscosity bound (1.8) is satisfied. Similar to a thermal state of 
M = 2* plasma with z/ = 0, we find that, while the bulk viscosity of the M = 2* plasma 
grows rapidly in the vicinity of the critical point for < z/ < 1 , it is finite precisely at 
T = Te(z/). 

We did not discuss in this paper gauge theory/string theory duality-motivated 
phcnomenological approaches to transport properties (see however [23]). We also did 
not discuss potential applications to RHIC/LHC physics (see however [24]). 

In the future it would be interesting to study transport properties of the Klebanov- 
Strassler [25] cascading plasma. The preliminary work necessary for such analysis 
already appeared in the literature [26, 21, 27]. 
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